Two-photon exchange and elastic electron-proton scattering 
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Two-photon exchange contributions to elastic electron-proton scattering cross sections are evalu- 
ated in a simple hadronic model including the finite size of the proton. The corrections are found 
to be small in magnitude, but with a strong angular dependence at fixed Q 2 . This is significant 
for the Rosenbluth technique for determining the ratio of the electric and magnetic form factors 
of the proton at high Q 2 , and partly reconciles the apparent discrepancy with the results of the 
polarization transfer technique. 
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The electromagnetic structure of the proton is re- 
flected in the Sachs electric (Ge(Q 2 )) and magnetic 
(Gm(Q 2 )) form factors. The ratio R = h p Ge/Gm, 
where /i p is the proton magnetic moment, has been de- 
termined using two experimental techniques. The Rosen- 
bluth, or longitudinal-transverse (LT), separation ex- 
tracts R 2 from the angular-dependence of the elastic 
electron-proton scattering cross section at fixed momen- 
tum transfer Q 2 . The results are consistent with R « 1 
for Q 2 < 6 GeV 2 Q, 0. However, recent polarization 
transfer experiments at Jefferson Lab [j| measure R from 
the ratio of the transverse to longitudinal polarizations 
of the recoiling proton, yielding the markedly different 
result R sa 1 — 0.135(Q 2 — 0.24) over the same range in 
Q 2 Q, which exhibits nonscaling behavior. In this letter 
we examine whether this discrepancy can be explained 
by a rcanalysis of the radiative corrections, in particular 
as they affect the LT separation analysis. 

Consider the elastic ep scattering process e(pi) + 
p(p2) e(j>3) +p(pi)- The Born amplitude for one pho- 
ton exchange is given by 

e 2 _ 

M = -i^u(p 3 h t ,u(pi) u(p 4 )T f "(q)u(p 2 ), (1) 

<r 

where the proton current operator is defined as 

r^) = F 1 (« 2 )y + i^^^%, (2) 

q — p4 — p2 = pi — P3 is the four-momentum transferred 
to the proton (Q 2 = —q 2 > 0), M is the proton mass, 
and Fi and F% are linear combinations of the Sachs form 
factors G E and G M (see Eqs. (fT5|) — (fTT?f) ) . 

The resulting cross section depends on two kine- 
matic variables, conventionally taken to be Q 2 (or r = 
Q 2 /AM 2 ) and either the scattering angle 8 or the virtual 
photon polarization e = (l + 2(1 + t) tan 2 (0/2)) . It 
can be put in the form 



da = A(tG 2 m (Q 



eG|(Q 2 )) 



(3) 



da = dao(l+8). Usually 8 is estimated by taking the one- 
loop virtual corrections of order a as well as the inelastic 
bremsstrahlung cross section for real photon emission. 

The LT separation technique extracts the ratio 
(Ge/Gm) 2 from the e-dependence of the cross section 
at fixed Q 2 . With increasing Q 2 the cross section is 
dominated by Gm, while the relative contribution of 
the Ge term is diminished. Hence understanding the 
e-dependence in the radiative correction 8 becomes in- 
creasingly important at high Q 2 . By contrast, the po- 
larization transfer technique involves a ratio of cross sec- 
tions, and is not expected to show the same sensitivity 
to the e-dependence of 8 0]. 

The amplitude Mi for the one-loop virtual corrections 
can be written as the sum of a "factorizable" term, pro- 
portional to the Born amplitude Mo, plus a remainder: 



Mi = f(Q 2 ,e)Mo+M! 
Hence to first order in a (a — e 2 /Att) 



<5 = 2/(Q 2 ,e) + 2 



Rej>4Mi} 
\M \ 2 ' 



(4) 



(5) 



where A depends on kinematic variables. This expression 
is modified by radiative corrections, expressed in the form 



The factorizable terms dominate, and include the elec- 
tron vertex correction, vacuum polarization, and the in- 
frared (IR) divergent parts of the proton vertex and two- 
photon exchange corrections. These terms are all essen- 
tially independent of hadronic structure. The hadronic 
model-dependent terms from the finite proton vertex and 
two-photon exchange corrections are expressed in M\. 
These terms are small, and are generally ignored |5(. 
The finite proton vertex correction was analyzed re- 
cently by Maximon and Tjon 0, who found 8 < 0.5% 
for Q 2 < 6 GeV 2 . It does not show a significant e- 
dependence, and so we drop it here. 

The factorizable terms can be classified further. Each 
of the functions /(Q 2 ,e) for the electron vertex, vacuum 
polarization, and proton vertex terms depend only on 
Q 2 , and therefore have no relevance for the LT separa- 
tion aside from an overall normalization factor. Hence of 
the factorizable terms, only the IR divergent two-photon 
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FIG. 1: Two- photon exchange box and crossed box diagrams. 

exchange contributes to the e-dependence of the virtual 
photon corrections. 

For the inelastic bremsstrahlung cross section, the am- 
plitude for real photon emission can also be written in 
the form of Eq. In the soft photon approximation 
the amplitude is completely factorizable. A significant 
e-dependence arises due to the frame-dependence of the 
angular distribution of the emitted photon. These cor- 
rections, together with external bremsstrahlung, contain 
the main e-dependence of the radiative corrections, and 
are accounted for in the experimental analyses Q{. 

In principle the two-photon exchange contribution to 
A4i, denoted .M 77 , includes all possible hadronic inter- 
mediate states (Fig.^). Here we consider only the elastic 
contribution to the full response function, and assume 
that the proton propagates as a Dirac particle. We also 
assume that the off-shell current operator is given by J5J , 
and use phenomenological form factors at the jp ver- 
tices. Clearly this creates a tautology, as the radiative 
corrections are also used to determine the experimental 
form factors. However, because S is a ratio, the model- 
dependence cancels somewhat, provided we use the same 
phenomenological form factors for both A4o and yVf 77 in 
Eq. ©. 

The sum of the two-photon exchange box and crossed 
box diagrams has the form 

Ar , _ c 4 f d*k lN a (k) N b (k)] f6) 
~ J (2tt)* [ Da (k) + D b (k)\ ' W 

where the numerators are the matrix elements 

N a (k) = u(p 3 )7m(^i _ ¥)lvu(pi) 

x u( Pi )T»(q - k)(tf 2 + # + M)T v {k)u(p 2 ), (7) 

N b (k) = u(p 3 )7„(j*3 + #)7/*u(pi) 

x u{ Pi )T»{q - k)W 2 + # + M)T v {k)u{p 2 ), (8) 

and the denominators are the products of the scalar prop- 
agators, 

D a (k) = [k 2 - A 2 ][(fc - qf - A 2 ] 

x[(p! - kf - m 2 }[( P2 + k) 2 - M 2 }, (9) 
D b {k) = D a {k)\ Pl ^ P3+k . (10) 

An infinitesimal photon mass A has been introduced in 
the photon propagator to regulate the IR divergences, 
and the electron mass m is ignored in the numerator. 



The implementation of Eq. 10 is the main result of this 
letter. However, we also want to compare with previous 
work, so a partial analysis of the leading terms in © is 
warranted. 

To proceed, we can separate out the IR divergent parts 
from the finite ones. There are two poles in the integrand 
of JSJl where the photons are soft: one at k — 0, and an- 
other at k = q. For the box diagram, the matrix clement 
can be written as the sum of a contribution at the pole 
k = plus a remainder, N a (k) = N a (0) + N a (k). Explic- 
itly, we have 

N a (0) = 4pi ■ P2 q 2 iMo/e 2 . (11) 

The matrix element at the pole k = q is the same, so 
N a (q) — N a (0) (this also follows from symmetry argu- 
ments) . This suggests that the dominant contribution to 
the box amplitude can be approximated as 

There are two assumptions implicit in this approxima- 
tion. The first is that the integral involving N a (k) is 
small, and contains no ultraviolet (UV) divergences from 
the F 2 part of the current operator (0) . Without hadronic 
form factors, Eq. J2J does in fact lead to UV divergences. 
We demonstrate below how to get around this difficulty 
by rewriting F\ and F 2 in terms of the Sachs form factors 
Ge and Gm- The second assumption is that the hadronic 
form factors have no significant effect on the loop inte- 
gral, and can be factored out. In essence, this assumes 
that the hadronic current operators occurring in Eq. © 
can be replaced by T A '(0) = 7 M for the vertex involving 
the soft photon, and by T^(q) for the other vertex. 

With these caveats in mind, the IR divergent box am- 
plitude from the pole terms can now be written as @ 

The four-point function arising from the loop integral 
has been evaluated analytically in the limit A 2 <C Q 2 
following 't Hooft and Veltman 0- 

A similar analysis of the crossed box amplitude shows 
that 

N b (0) = 4p 3 -p 2 q 2 iM /e 2 , (14) 

and hence 

In the lab frame (pi -p 2 = E\M and p 3 -p 2 = E^M), the 
total IR divergent two-photon exchange contribution to 
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FIG. 2: Difference between the model-independent IR diver- 
gent contributions of Eq. 1 (it and of the commonly used ex- 
pression H17H . 



FIG. 3: Difference between the full two-photon exchange cor- 
rection and the model-independent IR divergent result of Eq. 
GJ). 



the cross section is readily seen to be 



(16) 



a result given by Maximon and Tjon |6J. The logarith- 
mic terms in m cancel in the sum, while the logarithmic 
IR singularity in A is exactly cancelled by a correspond- 
ing term in the bremsstrahlung cross section involving 
the interference between real photon emission from the 
electron and from the proton. 

By contrast, in the standard treatment of Mo and 
Tsai (MT) the loop integral in (JT3J is approximated 
by setting the photon propagator not at a pole equal 
to l/q 2 . This results in a 3-point function K{— pi,P2) 
which, unfortunately, has no simple analytic form in 
the limit A 2 <C Q 2 ■ After a further approximation 
K{—pi,P2) ~ K(p 1: p 2 ), the total IR divergent result is 
given as y| 

S m (MT) = -2- (K( Pl , P2 ) - K(p 3 ,p 2 )) , (17) 

7T 

where K(pi,pj) = p % ■ pj dy ln(p 2 /A 2 )/p 2 and p y = 

Pty + Pj(i -y)- 

Because (Sir(MT) is the result generally used in ex- 
isting experimental analyses Q, 0, it is useful to com- 
pare the e-dependence with that of Sir. The difference 
<5ir — <5ir (MT) is independent of A, and is shown in Fig. [3 
as a function of e for Q 2 = 3 GeV 2 and Q 2 = 6 GeV 2 . 
The different treatments of the IR divergent terms al- 
ready have significance for the LT separation, resulting 
in roughly a 1% change in the cross section over the range 
of e. This effect alone gives a reduction of order 3% and 
7% in the ratio R for Q 2 = 3 GeV 2 and Q 2 = 6 GeV 2 , 
respectively. 

We return now to the implementation of the full ex- 
pression of Eq. ijfJJl. The full expression includes both 
finite and IR divergent terms (there is no need to treat 



them separately), and form factors at the 7p vertices. To 
avoid sensitivity to the UV divergences in the loop inte- 
grals arising from the F 2 part of the current operator J5J , 
we rewrite F\ and F 2 in terms of the Sachs form factors 



F 2 (q 2 ) = 



G E (q 2 )+TG M (q 2 ) 

1 + r 
G M {q 2 )-G E {q 2 ) 
1 + r 



(18) 
(19) 



Ge and Gm are taken to have the common form factor 
dependence Gsiq 2 ) — Gm(<? 2 )/Mp = G(q 2 ), with G(q 2 ) 
a simple monopole G(q 2 ) = —A 2 /(q 2 — A 2 ). We leave a 
fuller exploration of the hadronic model-dependence to 
a future paper. Effectively the F 2 part of the current 
then behaves like a dipole, and the loop integrals are UV 
finite for any choice of cutoff mass A. We have taken 
A = 0.84 GeV, consistent with the size of the nucleon, 
for which the results show a plateau of stability. The sen- 
sitivity to A is mild because the form factor dependence 
enters as a ratio in 5. 

The loop integrals in Eq. © can be evaluated ana- 
lytically in terms of four-point Passarino-Veltman func- 
tions |8j, and trace techniques used to implement the 
sum over Dirac spinors implicit in Eq. (0. This is a 
formidable task that is facilitated by the use of estab- 
lished algebraic manipulation routines. We used two in- 
dependent packages (FeynCalc and FormCalc [lojp. 
which gave identical numerical results. The Passarino- 
Veltman functions were evaluated numerically using the 
FF program [ll|. 

The model-independent IR divergent result of Eq. jTSjl 
is an appropriate benchmark with which to compare the 
full result 6f u i\. Because the IR behavior is the same, the 
difference <5f u n — Sib, is finite (i.e. independent of A). The 
results are shown in Fig. A significant e-dependence 
is observed, which increases slightly with Q 2 . The addi- 
tional correction is largest at backward angles (e — > 0), 
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FIG. 4: The ratio of form factors measured using LT sepa- 
ration (hollow squares), together with the global fit (dashed 
line). The unshifted LT data represent a binned average of 
all LT separated data points with normalization factors de- 
termined by the global fit in Ref. Q. Filled squares show the 
shift in the LT results due to the two-photon exchange cor- 
rections (offset for clarity) , and the solid line shows the effect 
on the global fit. Error bars have been left unchanged. The 
polarization transfer data 3j are shown as hollow circles. 



and essentially vanishes at forward angles (e — > 1). 

To consider the effect on the ratio R determined in 
the LT separation, we make a simplified analysis that 
assumes the modified cross section is still approximately 
linear in e. The results shown in Figs. |21 and [3] are com- 
bined, giving A = <5f u n — <5ir,(MT). For each value of Q 2 in 
the range 1-6 GeV 2 we fit the correction (1 + A) to a lin- 
ear function of e of the form a(l + be). The parameter b so 
determined behaves roughly like b w 0.0141n (Q 2 /0.65), 
with Q 2 in GeV 2 . For the LT separation, the corrected 
Eq. J3J becomes 



da = {clA)tG z m {Q z ) [l + {BR Z + 6)eJ , (20) 

where B = 1/(/i 2 t), and R is the corrected ratio R. Since 

a ~ 1, we have essentially R 2 = R 2 — b/B. 

The shift in R is shown in Fig. 0] together with the 
polarization transfer data. The effect of the additional 
terms is significant. Although some dependence on nu- 
cleon structure is expected, these calculations show that 
the two-photon corrections have the proper sign and mag- 
nitude to resolve a large part of the discrepancy between 
the two experimental techniques. Clearly there is room 
for additional contributions from inelastic nucleon excita- 
tion (e.g. the A + ). These have been examined previo usly 
in Refs. 01 m various approximations. Greenhut 0] 
used a fit to proton Compton scattering to calculate the 
resonant contribution to two-photon exchange, and found 
some degree of cancellation with the nonresonant terms 
at high energies. Further study of the inelastic region 



is required, including also the imaginary part of the re- 
sponse function |l3j . 

Direct experimental evidence for the contribution of 
the real part of two-photon exchange can be obtained by 
comparing e + p and e~p cross sections. (A4q changes sign 
under e~ — > e + , whereas M. 11 does not.) Hence we ex- 
pect to see an enhancement of the ratio a(e + p)/a(e~p) 
due to two-photon exchange (after the appropriate 
IR divergences are cancelled due to bremsstrahlung) . 
There are experimental constraints from data taken at 
SLAC 14] for E x = 4 GeV and E x = 10 GeV, which 
are consistent with our results. However, the SLAC data 
are from forward scattering angles, with e > 0.72, where 
we find the two-photon exchange contribution is < 1%. 
A more definitive test of the two-photon exchange mech- 
anism could be obtained at backward angles, where an 
enhancement of order a few percent is predicted. 
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